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Surface properties, temperature effects, cooling behavior. 


and observability of neutron stars have been studied. For this 
purpose the opacity of the surface layers is calculated both 
for a pure iron and a pure magnesium composition. It is found 


that the non-degenerate layers are only a few meters thick and 


in no case exceed 1% of the stellar radius. The star cools mainly 

8 o 

through neutrino emission when T ^ 1 ~ 4 x 10 K, but at lower 


temperatures the cooling is primarily through electromagnetic 
radiation. The neutrino cooling mechanisms included were the 


neutrino plasma process, the URCA process, and the neutrino 
bremsstrahlung process. The cooling behavior is quite compli- 
cated. The rate of cooling generally depends on mass, nuclear 
potential, and surface composition, among which the dependence 
on mass is the most significant. It will be hard to observe 


low mass neutron stars due to fast cooling rates. However, 





medium and high mass stars should still have temperatures 

exceeding about 2 x 10 6 °K on the surface for times of the 
3 5 

order of 10 to 10 years. Hence it should not be imposs- 
ible to observe massive neutron stars relatively close to 
us, if there is no x-ray emission of larger flux coming 
from the surrounding region. 

INTRODUCTION 

From studies of nucleosynthesis and stellar evolution 
it appears likely that the remnants of some supernova ex- 
plosions will probably contain a central core of highly 
condensed matter of the order of nuclear density whose main 
composition is neutrons, and of surrounding envelopes of 
ejected material which are expanding continuously (Arnett 1966, 
Cameron 1959 a, b, 1965a, Chiu 1964, Colgate and White 1965, 
Zwicky 1938, 1939, 1958). The central core, which may 
become as hot as ~10"^ °K at the peak, will cool down to 
~10^ °K within ~10^ sec or so (Chiu 1964) , due to the extremely 
high rates of neutrino cooling. Whether the 
remnant, if a stable neutron star, will ever be detectable 
depends critically on the cooling behavior of the neutron 
star itself and on physical conditions in the ejected en- 
velope . 

When the first, crude, information on the recently 
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discovered galactic x-ray sources became available (Bowyer, 

By ram, Chubb and Friedman 1964 a, Giacconi , Gursky, Paolini 
and Rossi 1962, 1963, Fisher and Meyerott 1964), it was 
tempting to identify these x-ray sources with neutron stars, 

g 

for a neutron star of photospheric temperature T e ~10 to 

7 o 6 

3 0 K would emit soft x-rays with peak near 30 ~ 3 A. if it 

can be regarded as a black body. The first series of cool- 
ing calculations showed that neutron stars of this range of 
temperature can be sufficiently luminous and last sufficiently 
long to be consistant with the early observational data (Chiu 
and Salpeter 1964, Morton 1964, Tsuruta 1964). In the mean- 
time, possible cooling mechanisms have been reexamined and 
some faster cooling mechanisms have been proposed (Bahcall 
and Wolf 1965 a, c, Finzi 1965a, Tsuruta and Cameron 1965 a). 

The results of our latest calculations (whose details are 
given in this paper) indicate that a neutron star with 

T ~10 7 °K will cool too fast, 1 day to 10 years (for models 
e 

of different mass) , to be detected as x-ray sources, but that 
6 o 

if T ~2xl0 K, the medium and high mass stars will last 
e 

3~5 

sufficiently long for detection (~10 years) . The conclusion 
is that the cooling rate alone does not exclude the possibility 
that some neutron stars will be detectable, provided that the 
thermal soft x-ray emission exceeds the nonthermal background 
from the ejected envelope. 
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If we look at a bare neutron star in a vacuum, there 
would be several characteristic features which ecu Id be checked. 
For instance, a neutron star (of radius ~ 10 km) should appear 
to be a point source, and the x-ray spectrum from it should 
approximate that of blackbody radiation. A lunar occultation 
measurement conducted by Boyer, Byram, Chubb and Friedman (1964 b) 
showed that the source in the Crab Nebula has a diameter of 
about one light year. Clark (1965) and others reported that 
the same source emits a significant x-ray flux in ~ 30 kev 
region, which is inconsistant with the spectrum of a blackbody 
radiation from a surface of a few million degrees. The MIT 
and American Science and Engineering group and the Livermore 
group reported that the spectrum of the strongest x-ray source 
in Scorpius is inconsistent with the hypothesis of blackbody 
radiation (Giacconi, Gursky and Waters 1965, Chodil, Jopson, 

Mark, Seward and Swift 1965) . These observational results 
seem to indicate that the observed x-rays from some sources 
do not consist alone of blackbody radiation from a neutron 
star. However, we cannot conclude that these results are 
evidence against the existence of neutron stars. 

One of us (Cameron 1965 b) suggested that a possible 
model of an x-ray source may be a vibrating neutron star 
with an associated magnetosphere and surrounding hot 
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gas clouds. (Also, see Finzi 1965 b) . The central 
neutron core may emit thermal x-rays and the 
surrounding area may emit non-thermal x-rays (either 
synchrotron radiation or bremsstrahlung) . The question of 
whether the observed x-rays are thermal, non-thermal, or 
a mixture will require more extensive observations. Some 
of the recent experiments quoted above only suggest that 
the thermal component is not the main contribution to the 
observed x-rays, and that the effects from the surrounding 
envelopes are dominant. The problem of whether we can dis- 
tinguish between the thermal component from the neutron 
star and the non-thermal components from the surrounding 
area will be considered in the last section of this paper. 

As long as there is a possibility that neutron stars are 
related to the galactic x-ray source, directly or indirectly, 
the study of neutron stars will continue to be of astronomical 
importance, as well as of fundamental importance in general 
relativistic physics. 

In our earlier paper (Tsuruta and Cameron 196 6 ) pro- 
perties of neutron stars which are independent of temperature 
were studied. Here we discuss some thermal properties of 
these neutron stars. For this purpose six typical models 
have been chosen from our first paper. Three of these are 
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based on the V p type of nuclear potential, and the other 

three are based on the V type of potential. The different 

v 

types of potential correspond to different assumptions about 

the possible nuclear forces which will govern the interactions 

between the constituent baryons . Among the three models for 

each type of potential, one is of low mass at the low density 

base of the principal mass peak, one is half-way up the peak, 

and the other is near the top of the peak. The characteristic 

properties of these models are listed in Table 1. Each model 

is designated by its potential type and approximate mass value. 

It has been suggested that a model with density p > 8 P nuc i^P nuc l 

14 3 

is the nuclear density=3.7 X 10 gm/cm ) is unreliable (Bahcall 
and Wolf 1965 a, b) due to the uncertainty in high energy 
physics. Our model (V^,1M 0 ) is denser than this critical limit. 
However, all other models lie roughly within the limit of val- 
idity and the general conclusions deduced from these models are 
expected to be reasonably reliable. 

Typical models of quite different mass value and nuclear 
potential have been chosen in order to avoid as far as possible 
the danger of drawing general conclusions from limited assump- 
tions. The main approach was, therefore, to set upper and 
lower limits which define the range where the most probable 


models of neutron stars would lie. 



In the following sections the surface properties of 
neutron stars will be studied first, the cooling and re- 
lated calculations will be presented, and finally the 
problem of detectability will be discussed. 


ENVELOPE STRUCTURE EQUATIONS 
The structure equations for the outer layers of a 


neutron star are : 
dP 
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For convective equilibrium: 
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where M , 
r 

T r' P r' p(r) and 

k (r) 

are the mass, temperature. 


pressure, density and opacity at a distance r from the center, 
L is the total luminosity of the star, G is the constant of 
gravitation, a is Stefan's radiation constant, c is the vel- 
ocity of light, and F is the ratio of the specific heats 
c^/c^. The assumption was made that the star is spherically 
symmetric and in hydrostatic equilibrium and that there is no 


energy generation or sinks in the envelope. The general 
relativistic expressions are used in the hydrostatic equa- 
tions (1) and (2) because the general relativistic effects 
are quite important even near the surface for some of the 
denser stars (Tsuruta and Cameron 1966). 

Near the surface, the above equations are more easily 
integrated in electronic computers if they are expressed 
in logarithmic form, and it is better to use pressure as 
an independent variable because the pressure gradient is 
quite high near the photosphere of neutron stars . Then 
Eqs. ( 1 ) — ( 3 ) can be expressed as: 

d ( ^n r ) _ exp ('knP 1 -tr A- -In b- £n g- tn d) 
d(^np r )“ ~ r 

d(^nMr) ^ ^ exp ( 3 inr-im x +ln(4£yinp (r) 
d ( £ nP G 
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r 

d (^ nT r) _ exp [tn(3/16TTacG) +ln (r)-!-hi P (r) hriL-4-t rff -hi 
d (tnP r 

(7) +hnP +hnA-tnB-tnD] 

r 

where 

A=exp (hnr)-exp {tr\2+^rM -2hnc+tnG) 

r 

(8) B=exp (tn.p (r J+exp (: p nP r -2hnc) 


(5) 


( 6 ) 


D=exp (-InM^J+exp (tn4ir+3t nr+^nP r -2^nc) 
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The equation of state in non-degenerate layers is 


( 9 ) 


=(i 1 ^ 1 4 

v— + — y — + — a T 


^ion 


+ — a T 
H 3 


where the first term is the gas pressure and the last term 
is the radiation pressure. In degenerate envelopes, all 
pressures except the degenerate pressure of electrons are 
negligible and the following equation for degenerate electron 
gases is applicable: 


p *n.U e H = B X , 


( 10 ) 


P=A f (x) ; 
f (X)=X(2x"-3) (x"+ 1) 2 +3 sinh _J -X 


X=P F /m e c 


2 _ . 2 , _ . ,-l 
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A=nm c /(3h ) =6.01x10 , B=8nm c q H/3h =9.82x10 q . 

e e e e 


d and u. are defined as 
e ion 


p = [Z A n(A Z )]/ il Z n (A Z )] 

e jfa -vt it / v i 1 i/i 
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ion 


= [E A.n(A.,Z») ]/[E n(A. ,Z. )] 
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where n(A^#Z^) is the number density of the nucleus i of 
mass number A^ and atomic number Z^ , H is the proton mass, 
and the remaining notation is conventional. The equilibrium 
composition of matter as found by Tsuruta and Cameron (1965b) 
was used. In the envelopes of neutron stars the equation of 
state is relatively simple as shown above, but the behavior 
of the opacity is quite complicated. For this reason, the 
opacity is treated separately in the next section. Generally, 
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great care has to be taken in evaluating P if the envelopes 
are in convective equilibrium. However, as is shown 
later, convection appears to play no role in neutron 

stars . 

To obtain a proper boundary condition at the surface, 
it was assumed that the ordinary theory of stellar atmospheres 
would apply in the atmospheric layers above the surface of 
neutron stars, provided that general relativity effects are 
correctly taken into account in some of the denser models. 

The surface of the star is defined as the point where the 
actual temperature is equal to the effective temperature T^ 
(the temperature of the black body which would radiate the 
same flux as the star itself) . Then the total luminosity of 
the star L is expressed as 

2 4 

(12) L = 4ttctR T 

' ' e 

a is Stefan's constant. 

If we assume that the opacity is independent of both height 
and wave length in the atmosphere but that it has the constant 
value determined at the photosphere, the theory of radiative 
transfer in stellar atmospheres leads us to the following 
simple relation: 


(13) 
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where 


is the optical depth at the surface, (p 
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opacity at the photosphere and the general relativistic form 

“1 


of g(p v,, p„u) is: 

ph ph 2 p 

GM "P i1l ,Y, 4nR P h y, 2GM 

(14) 9= — 2 ! _Eh_ 3+ 1,^1+ 5 X, 1 - 2 

K p ph c m rc 


). 


The subscript ph stands for the value of the respective vari- 
able at the photosphere. Noting that the mass content and the 
thickness of the atmosphere of neutron stars are negligible, 
the radius and the gravitational mass of the cold models of 
neutron stars can be used for R and M above. Then equations 
(9) , (13) and (14) evaluated at the photosphere give us 

sufficient boundary conditions at the surface of a neutron 
star of surface temperature T . 


OPACITY 

The total opacity k in a neutron star can be expressed as 

11 . 1 

(15) K + K 

R c 

where k and h are the radiative and conductive opacity, 
respectively. Radiative opacity is due to the various pro- 
cesses of atomic and molecular absorption, emission, and 
scattering of radiation in which electrons play the major 
role. The relative importance of these processes depends 
strongly on the temperature-density combination. For in- 
stance , in matter of high temperature and of relatively low 
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density, electron scattering is dominant, while in the region 
of intermediate density and temperature the various photo- 
electric effects are the most important. In degenerate matter 
of high density, electron conduction is the most efficient 
mechanism. The major processes of atomic absorption are 
the bound- free, free- free, and bound-bound processes. Excited 
electrons emit photons in the inverse processes. The scatter- 
ing processes are Thomson scattering if T<5xl0^ °K and Compton 
scattering for higher temperatures. 

The opacity generally depends on density and temperature 
in quite a complicated way. In recent years, various ex- 
tensive tables based on detailed computations have been 
published which give the absorption coefficient for many 
different compositions and for a large number of points in 
the temperature-density diagram. The most accurate method 
of obtaining opacities at present appears to be through use 
of the computer code for opacities constructed by A.N.Cox 
and his colleagues of the Los Alamos Scientific Laboratory 
(Cox, 1961) , which includes most of the possible major 
processes contributing to opacity, and this code was kindly 
made available for our calculations. It includes bound- 
free and free- free abosrption, electron scattering (both 
Thomson and Compton scattering) , negative ion absorption 
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and electron conduction. The bound-free absorption depends 
on the equilibrium number of electrons which are bound in 
the various atomic states . When the ionization of one ele- 
ment is completed, no more bound- free absorption due to that 
element can occur. For high densities the effect of de- 
generacy is taken into account in all but the electron 
scattering term. At low densities and low temperatures not 
all electrons are ionized. An ionization code was used in 
conjunction with the opacity code in these regions to cal- 
culate the degree of ionization, the partial pressure of 
electrons, and the number of free electrons in the opacity 
calculations . 

At temperatures above about 5xl0 7 °K, almost all ele- 
ments are ionized. The existing tables are used to obtain 
absorption cross sections for various kinds of processes. 

The electron scattering term is obtained for the non-degen- 
erate case, and pair production of electrons and positrons 
is not considered. Therefore, the opacity is independent 
of density but dependent on temperature in the high tempera- 
ture region where Compton scattering is dominant. Different 
approximations are applied for different degrees of de- 
generacy to evaluate the conductive opacity. 

Values of the opacity were calculated for a pure iron 
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composition and a pure magnesium composition, in the range 

3.7 o 10 o 

of temperature from 10 K to 10 K, and of density from 

-43 3 14 3 

10 * gm/cm to 10 gm/cm . The reason for these particular 

choices of composition is explained in the next section. 

14 3 

Opacities at densities higher than 10 gm/cm have not been 
included because degeneracy sets in at densities far below 
this. Calculations at temperatures higher than 10^ °K has 
not been carried out because the assumptions of nondegenerate 
electron scattering and no electron-positron pair creation 
cause serious errors there. Also, neutron stars of tempera- 
ture higher than this are of no interest to us because they 
would cool too quickly. The case T^ <10"^*^ °K and p <10 ^ ^ 
gm/cm" 3 was not included because the opacity code did not 
work in these low temperature, low density regions. 

Results obtained are plotted in Figure 1. The solid 
curves represent the opacity of iron 56 as a function of 
density at different temperatures, while the dashed curves 
represent the same for magnesium 24. The opacity shows quite 
a complicated dependence on density and temperature in the 
region p s^lO 6 gm/cm" 3 and T<10 8 °K, where the transition from 
electron scattering or electron conduction to bound- free 

opacity occurs. The almost straight horizontal lines for 
9 o 6 3 

T^10 K and p^lO gm/cm are due to Compton scattering. 

The almost straight lines of negative slope in the region 
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of high density are an indication that electron conduction 
is the dominant factor in the transport of energy there. 

We can assume that degeneracy starts as soon as the opacity 
in Figure 1 starts to follow one of these straight negatively 
sloping lines. The opacities (Inn) thus obtained have been 
stored as an input deck of cards in the form of a two-dim- 
entional table corresponding to given InT and lnp combina- 
tions, for later use. 

The results from the opacity code calculations were 
checked in the various asymptotic regions of density and 
temperature, using simpler, analytic approximations. In the 
region where photoelectric effects are dominant, the Kramers 
opacity formulae (see, for instance, Schwarz schild 1958) 
were used. The opacity due to electron scattering was checked 
through the equations given by Sampson (1959) . In order to 
check the opacity in the region of electron conduction the 
relations in Schatzman (1958) were used. The agreement was 
quite satisfactory. It turned out that the neglect of degen- 
eracy and electron-positron pair creation in the formulae 
for electron scattering in Cox's opacity code causes no 
serious errors in the problem of neutron stars. 
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ENVELOPES AND THEIR CHARACTERISTICS 

A computing program has been prepared for the 7094 com- 
puter which carries out the surface integrations and other 
related computations automatically. The program was constructed 
so that the equation of state automatically switched over from 
Equation (9) to (10) as soon as the point was reached where 
these two became equal. The surface boundary values were 
calculated through the subroutine for the atmosphere whenever 
a new set of values of radius R, mass M and photospheric 
temperature T of the star were given. The interval of in- 
tegration Ain P was automatically adjusted so that the change 
of every variable was kept smaller than a suitable preassigned 
limit. For a given In T and In p combination the correspond- 
ing opacity (Inn) was obtained by linear interpolation in the 
input opacity table. The surface boundary values of two 
typical models of stable neutron stars (V^,0.6 M q ) and(V^,2 

with a pure iron atmosphere are shown in Table 2 . 

To determine the temperature distribution in a typical 

neutron star envelope, the integration was first carried out 

14 3 

from the photosphere down to the point where p=10 gm/cm , 
for a representative model of M=1 and R=10 km. This was 
repeated at several different surface temperatures. The 
results are shown in Table 3. The temperatures at different 
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densities are listed in terms of the given surface tempera- 
ture T^. and are the temperature and density where 

degeneracy starts (where equation (9) gives the same pressure 

6 3 

as equation (10)). Degeneracy starts at about p^=10 gm/cm 

when T ~10 7 °K. But when the surface has cooled down to 
e 

about 10 6 °K, degeneracy sets in at p^=10^ gm/cm ^ . A signi- 
ficant result is that even after the degeneracy boundary 
has been passed, the temperature still continues to go up as 

we go inwards. The fractional rise in temperature as the 

6 9 3 

density increases from 10 to 10 gm/cm is about 10% when 

T e =7.7xl0~* °K, but at T^IO^ °k the temperature at the point 

p=10^ gm/cm^ is about 3 times that at p=10^ gm/cm^. As we 

9 3 

go in toward the center from the point with p=10 gm/cra to 
12 3 . 

p=10 gm/cm , the fractional rise m temperature is about 

0.5% for T =10^ °K and about 3% for T =10 7 °K. The table 
e e 

shows that the temperature gradient is completely negligible 

12 3 

for density higher than about 10 gm/cm . The conclusion 
is that the temperature gradient is very high in the outer- 
most thin non-degenerate layers, the temperature continues 
to rise as we go through the degenerate outer layers of 

heavy ions and electrons, but the inner neutron core (with 
12 3 

p^ 10 gm/cm ) is isothermal even for the models of the 


hottest neutron stars. Hence, the core temperature T (which 

c 
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is also the central temperature of the star) can be defined 

12 3 

as that temperature where p=10 gm/cm . 

The central temperatures are plotted against surface 
temperatures in Figure 2 . The solid curves are for Fe and 
the dashed curves are for Mg atmospheres. Curves drawn for 
three models of the type are marked by the corresponding 
mass. Similar but simpler calculations were made by Chiu 
and Salpeter (1964) and Morton (1964) . Their central tempera- 
ture was defined to be the temperature where the degeneracy 
sets in (our T ) , but we have noted that the temperature 
continues to rise considerably as we go inwards passing the 

degeneracy boundary (T >T, ). Therefore, our central tempera- 

c b 

ture should be higher than their values in general. However, 

in high temperature regions where the electron scattering is 

the most important mechanism for the opacity, they used the 

2 

constant value h~ 0.2 cm. /gm, Thomson scattering opacity, 

while the Compton scattering included in Cox's code which 

we used lowers the opacity from the constant value of Thomson 

scattering. This, in effect, lowers our values of central 

temperature. These two causes of discrepancy compensate one 

another and there is good general agreement among our results 

9 o 

and their results for T c >10 K. We see in Figure 2 that the 
central temperature is somewhat lower for Mg than Fe at the 
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same surface temperature. This is due to the fact that some- 
what lower opacities are associated with Mg than Fe as is 

revealed in Figure 1. For some of the coolest neutron stars 
4 o 

(T^~10 K) the central temperature is only about 10 times 

the surface temperature, while for hot models ( of T ~10 7 °K) 

e 

the core is about 100 times as hot as the surface. In any 
case, however, the ratio of the central temperature to the 
surface temperature is quite small as compared with that of 
ordinary stars. 

To examine the region near the surface more in detail, 
temperature is plotted against distance from the photosphere 
as measured inward in Figure 3. Each curve is marked by the 
surface temperature. The model with one solar mass and 10 km 
radius is used to illustrate the general behavior of the 
surface properties. The crosses marked by y=2.5 represents 
points where the degeneracy starts. This criterion for de- 
generacy is derived from the fact that the kinetic energy of 
a non-relativistic fermion (about 3/5 of the Fermi energy E ) 
and the thermal energy of a free particle with no internal 
degrees of freedom (3kT/2) should be equal at the boundary 
between the non-degenerate and degenerate layers. The re- 
sult of the present calculations shows that even for the 
hottest models degeneracy starts before we go inward by 100 
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meters from the surface and that the non-degenerate layers 

are less than 1 % in thickness for even the hottest models. 

The mass contained in the non-degenerate envelopes is 

practically zero. In our previous paper (Tsuruta and 

Cameron 1966) it is seen that the amount of mass contained 

even in the inner degenerate electron- ion envelopes is very 

small compared to the total stellar mass. These results 

more than justify our previous assumption of constant mass 

and radius in the atmospheric calculations and also the 

neglect of non-degenerate layers in determining the total 

mass and radius of the star in our previous paper (Tsuruta 

and Cameron 1966) . Hot neutron stars with T ~10^ °k have 

e 

non-degenerate envelopes of about 10~20 meters thick but 
when the surface temperature falls to about a million degrees 
the thickness of the non-degenerate layers becomes only about 
a meter or so. A typical neutron star with the sun's mass, 

10 km radius, and 6.7x10^ °K at the surface (with ~100L ) is 

G 

shown to have non-degenerate envelopes of 3~4 meters. 

The density profile near the surface is plotted in 
Figure 4 for the same model. The distance from the surface 
is now shown in centimeters. Within about a meter (0.01% of 

the radius) from the photosphere, the density rises to about 

5 3 6 o 2 5 3 

10 gm/cm when T e ~10 K and to about 10 * gm/cm for 
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hotter stars (T -1.6x10 K) . In the photosphere the density 
rises within a thickness of 10 cm by a factor of about 100 

^ Q 

for typical models (T ~10 K) . In any case a sharp drop 

14 15 3 

of density from the central value (10 —10 gm/cm ) to 

3 

the photospheric value (0.01—1 gm/cm ) occurs only at the 
very edge of the star . 

The distribution of density, temperature and degree of 
degeneracy E /kT within the thin layers about 10 meters from 
the surface are numerically shown in Table 4 at several in- 
teresting values of surface temperature. On comparing this 
table with the previous one, we see that the degeneracy 
criterion used in these two tables agrees well with each 
other. 

The conclusion according to the present calculation 

3 

is that neutron stars of about 10 times solar luminosity 

are as hot as 10 °K at the surface and about 10^~‘*'^°K in the 

interior, those which are as luminous as the sun (in the 

6 o 

x-ray region) are about 1—2x10 K at the surface and about 
8—9 

10 °K in the interior, and that by the time they cool down 
to the point where T^~ 10^ 7 °K and T^— 10^ °K they are too 
faint to be detected (L— 10 ^ Lq) . In our previous papers 
(Tsuruta and Cameron 1965b, 1966) it was shown that the composi- 
tion of the surface layers changes sharply from layer to 
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layer. Starting from the boundary between the neutron core 

and the degenerate electron- ion envelopes the composition 

changes from more neutron-rich heavy nuclei to less neutron- 

rich ones as we go outwards. In the outermost non-degener- 

6 3 

ate layers with P<10 gm/cm , the main equilibrium composi- 
tion should be ordinary iron group nuclei. This is why 
iron was chosen in our opacity calculations. 

A possible change of surface composition can occur if 
a diffusion process is present. The diffusion process can 
become quite efficient in the presence of very small density 
scale heights and large gravity effects, as is the case in 
the atmospheres of neutron stars. Rough estimates of the 
effect of diffusion on the surface composition of neutron 
stars were made by Chiu and Salpeter (1964) . Their con- 
clusion is that some lighter elements such as Mg, 0 and Ne 
can be present on the surface of neutron stars. This is 
why not only Fe but also Mg was selected in our opacity 
and atmospheric calculations earlier. The change of com- 
position will not occur if convective mass motions in non- 
degenerate layers cause efficient mixing of elements. In 
this case the original statistical equilibrium composition 
of iron will be maintained. However, convection appears 
to play no important role in neutron star problems since 
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the temperature gradients in our model atmospheres are all 
subadiabatic . At the present time, the effect of the dif- 
fusion is not known. However, the difference between our 
results for Fe and Mg compositions is relatively small. 
This indicates that the uncertainty of surface composition 
due to the effect of diffusion will not cause any serious 
errors in our results presented in this paper. We will 
see that the uncertainties due to other effects are far 
greater . 

ENERGY CONTENT OF A NEUTRON STAR 


If we assume that a neutron star belongs to the end 
state of a thermonuclear evolution, there can be no energy 
generation within it. Any stable neutron star is already 
so dense that gravitational potential energy due to con- 
traction is not available. Even though the matter is 
highly degenerate, the only contribution to the total 
energy of a stationary neutron star comes from the small 
tail of the Fermi distribution function of the particles 
which constitute the star. The heat capacity of a nearly 
zero-temperature ideal Fermi gas was derived by Chandrasekhar 
(1957). Using his result, the total thermal energy of a 
neutron star was found to be 


2 7 

r /'tt ^ 

u=l ’oi ^ 2 


(X^+l) 1 / 2 ^ 2 

5 — y 2 1 )n. 4rrr dr 

\ y i 
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F 

where X-; = P.:/^ 0 ) 

1 l 

F 

m i / ^*i an< ^ n i arS tlle mass ' Ferm;L momentum and number density 
til 

of the i particle, R is the stellar radius and T is the 

temperature of the isothermal core. The summation was taken 

over all the baryons and leptons present in each layer dr. 

The total energy was calculated as a function of temperature 

for our models. The results are shown in Tables 5-10. It 

is seen that as the surface temperature of a star decreases 

from about 5 x 10^ °K to 10^ °K the energy content of the 

50 40 

star decreases from about 10 ergs to 10 ergs , although 
the precise value depends on the type of model in question. 

In the above derivation the nuclear interaction between 
particles was neglected. The presence of nuclear forces 
will modify the heat capacity by an amount which typically 
can be of the order of factor 2 , according to rough estimates 
which we have made. We did not take such modifications into 
account in making the actual cooling calculations, because 
the uncertainties arising from other sources are far greater. 

NEUTRINO LUMINOSITY 

The conserved vector current theory of weak interactions 
predicts that neutrino- antineutrino pairs can be radiated 
in quantum electrodynamic processes as well as electromagnetic 
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radiation. Even though the probability for the neutrino 
radiation is enormously small, it plays an extremely import- 
ant role in some stages of stellar evolution because of the 
fact that the neutrino mean free path is so large that it 
can escape even from a dense star with hardly any inter- 
action, while electromagnetic radiation can only diffuse 
out very slowly from the interior to the surface. Various 
different neutrino processes possible in a stellar interior 
have been examined. Consequently, the following three 
processes have been found to be the most important in the 
problem of neutron stars. 

(1) Neutrino pair emission from the plasma process: 

(17) v (plasmon) -*v +v 

e e 

These neutrinos arise from the decay of plasmons in the 
degenerate electron gas in the interior of the neutron star. 
The rates are given by (Adams , Ruderman and Woo 1963 , Inman 
and Ruderman 1964) ; 

3 22 /'flu) *\9 

Q (ergs/cm -sec) =1 . 228x10 ( — °~2 ) F(x) 

pi \m e c 

where F(x)=T , K 2 or 

n=1 “nx 

X 3 F (X) =2C (3 ) +ix 2 £nx-i (2£n2+l) X 2 +X 4 fny 

96 


(^n 2-i+fn 2tt - '^(j) 0 ^ for X<2tt 


C(3)=| =1 (^;C(2)=| =i e i Vf- ; - C'(2)=| =1 

n n 


ri \ n 


fnn 


n 


,., 2 -„„2 r ,A> j'l-lM ~ 4 e 2p F 3/( 3 „ E ^) 


X ^t' < =4re "J d ' P E 


t 

E 
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E f - < P F + ^ ■ P f ■ < 3 " 2 n e> 1/3 

for relativistic degenerate electrons. (In the last three 

equations the units are = c = 1.) Then the neutrino 

luminosity due to the plasma process is 

( 19 ) = r R 0 4nr 2 dr (ergs/sec) 

v J o pi 

To evaluate the above integral, a table listing F(\) as a 
function of x and the interpolation subroutine were used. 
The plasma neutrino emission rates have been calculated as 
a function of temperature and electron number density and 
are shown in Figure 5. 

(2) The URCA process (the beta process and its inverse 
process) s In the interior of a neutron star it can be 
represented by reactions such as: 


(20a) 


n + n- n + p + e +v 


(20b) 


n + n-»n + p + fi + v 


(20c) 


n + p + e -.n+n + v 


(20d) 


n+p + p, -»n + n + v 


The approximate formulas for the rates of these processes 
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for a neutron star of uniform density distribution have been 
derived by Bahcall and Wolf (1965c). We applied their 
results to our neutron star models „ The neutrino luminosity 
due to the URCA process can then be given by 

L = j 10 (l+F)T g (p/P nucl ) 4 nr dr (ergs/sec) 

(21) F = P F (n)/ P p (e) = [1-2.25 (P nucl /0) 4/3 ] V2 for p>1.8p nuel 

=. 0 for p < 1.8 p , 

nucl 

where p is the density of the neutron star matter, p , is 

nucl 

14 3 

the density of nuclear matter (= 3.7 x 10 gm/cm ), T Q 

is the stellar temperature in units of 10 °K,P i: ,(i) is the 

th 

Fermi momentum of the i species. 

(3) The neutrino bremsstrahlung process: In the interior 

of a neutron star it can be symbolically expressed as 

mm -J- mm -J- 

(22) e + (baryon)— - e + (baryon)— + v + v 

G G 

The neutrino pairs are emitted when electrons scatter from 
positive or negative baryons in the interior of the neutron 
star. Ruderman and Festa (private communication) have 
kindly provided us with the following approximate preliminary 
expression for this process: 
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(23) Q (ergs/gm-sec) = 10 6 Z 2 — (T ) 6 for E » me 2 
H n y r 

where Z is the effective charge of the electron scattering 

centers, n is the number density of such centers, and n 
z 

we take here to be the baryon number density. Ruderman and 

Festa have suggested that there may be proton clustering in 

neutron star interiors in the presence of very large numbers 

of neutrons, so that the effective charge of a scattering 

center might be 2. However, for some of the massive neutron 

15 3 

stars in which the density exceeds about 10 gm/cm there 
will also be a large number of Z hyperons, which may hinder 
the clustering process., In such a case, we have chosen to 
take the effective charge of a scattering center to be 1 and 
to count as the scattering centers both the protons and the 
T~ hyperons. The possible error caused by the uncertainty 
due to the different interpretation of the scattering centers 
and the effective charge is in any case negligible in our 
present problem* The neutrino luminosity due to the brems- 
strahlung process is 

B r R 2 

(24) L = ; Q p 4rr r dr (ergs/sec) . 

V < B 
o 

The total luminosity is obtained by adding all the 


contributions 
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L = L . + lP l + l URCA + l B 


ph 


V 


V 


is the photon luminosity defined by (12). 

The results are shown in Tables 5-10. To visualize the 

contributions of different processes, the three different 

neutrino luminosities and photon luminosity, together with 

the total energy, are plotted against the central temperature 

for the model (V , 1.1M) and shown in Figure 6. At higher 

Y © 

temperatures the URCA process and the plasma neutrino process 

compete with each other but the URCA process is always the 

most important „ At lower temperatures the URCA process 

competes with the bremsstrahlung process but the plasma process 

becomes unimportant. At sufficiently low temperatures the 

bremsstrahlung process predominates over other neutrino 

processes. The point where the bremsstrahlung rate begins 

to exceed the URCA rates depends on the stellar mass. For 

g 

massive stars, this switching occurs at 3^8x10 °K, but for 

7 

low mass stars the switching temperature is ~ 5 x 10 °k„ 
Neutrino cooling of any kind becomes too small as compared 
with the photon cooling when the temperature becomes lower 

3 

than about 1~ 4 x 10 °k. Hence, the branssfcrahlung process 
never becomes important in sufficiently low mass stars 

(M ~ 0o 2M ) . 

o 
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In the above luminosity calculations we assumed that all 
neutrino processes not included here are unimportant in 
neutron stars. Bahcall and Wolf (1965a, c) have raised the 
question of neutrino cooling from pion decays in neutron star 
interiors. Such pion decays can occur only if pions should 
have a small effective mass in the presence of a largely 
neutron gas. Both Bahcall and Ruderman, as well as others, 
have indicated to us (private communication) their expectation 
that, under the conditions in which pions may be present in 
a neutron star, there will be a predominantly repulsive 
interaction between the pions and the neutrons. This would 
raise, rather than lower, the effective mass of the pions, 
and make it very unlikely that pions will be present in the 
interiors of neutron stars. 

COOLING TIMES 

The cooling time t is computed from the various theore- 
tical data of the last sections and from the following 
relation : 


r 2 j™. 

jy L(U) 


(26 ) 
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where U is the total energy content of a star and L (U) is 

the total luminosity, L + £ L , expressed as a function of 

ph v 

U 0 If the above integration is carried out from the initial 
energy to the final energy t gives the time interval 

during which the star has cooled from a higher temperature 
where the total energy is to a lower temperature where 
the total energy is U^. The moment at which a supernova 
explodes has been defined as the starting point for counting 
the age of a neutron star. 

The results are tabulated in Tables 5 through 10. An 
important parameter which characterizes the cooling behaviour 
is the ratio of the central temperature to the surface temp- 
erature, a typical value of which is about 100 for a neutron 
star. Hence we defined a as the ratio of the central temp- 
erature to the surface temperature in units of 100, and the 
values of a for each model at different temperatures are 
listed in the same tables. We see that this ratio is a 
sensitive function of mass, temperature, composition and 
nuclear potential. With increase in mass, a decreases. With 
increase in temperature, a increases. For the models of the 
same mass, temperature and nuclear potential, the value of 
a for Fe models is somewhat larger than that for Mg models. 
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We see that the detailed behavior of cooling is different for 
each different type of model and that this difference comes 
mainly through the sensitive dependence of a on the different 
kind of model. 

Figure 7 shows the cooling curves of our six models with 
iron atmospheres (solid curves correspond to type and 
dashed curves to type models). The point at which the 
neutrino cooling rate and the photon cooling rate become equal 
is shown by a cross for each curve. The figure shows that for 
the same nuclear potential, the neutron stars of lower mass 
cool faster than those of larger mass up to an age of about 
10 years, but after that the heavier neutron stars cool 
somewhat faster. 

The complicated effect of the nuclear potential is 
observed when we compare the curves of the models of the 
same temperature and mass but of different nuclear potential. 
To see better the effects of different composition, cooling 
curves for the same types of model of the same masses but of 
different compositions are shown in Figure 8. Solid curves 
represent models with iron atmospheres and dashed curves 
those with Mg atmospheres. We see that in the region where 
the neutrino cooling predominates over the photon cooling 
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the cooling rates of Fe models are somewhat faster than the 
cooling rates of Mg models, but that the reverse situation 
is noted in the region where the photon cooling is the main 
cooling mechanism. This is easily explained if we note that 
the neutrino cooling and the total energy content of the star 
depend on the internal temperature while the photon cooling 
is a function of the surface temperature and that the opacity 
of Mg atmospheres is somewhat lower than the opacity of Fe 
atmospheres c 

7 

We see in these figures that the models with T^ > 10 °K 

cool too fast for observation, while the strong absorption of 
x-rays from a neutron star by interstellar gases makes it 
very difficult for us to observe neutron stars with < 10^ °K<, 

rr 

The star's luminosity itself is already low at T g ~ 10° °K 

(Table 2) „ Hence the most important range of temperature of 

neutron stars from the point of view of observation is 
7 6 

10 > T > 10 °K, and this portion of the curves is enlarged 

e 

in Figure 9o A neutron star will be only about 1 day to 10 

*7 

years old when T g ~ 10 °K depending on its mass value, 

nuclear potential and surface composition, but it will take 
3 5 

about 2 x 10 to 3 x 10 years before it will cool down to 
£ 

~10 °K. We note that the effect of mass on cooling time is 
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the most drastic. The dependence of cooling on the different 
possible kinds of nuclear potential and composition is by no 
means negligible when we need detailed theoretical informa- 
tion, but this uncertainty is relatively small as compared 
with the mass effect. By comparing our present results with 
the original calculations (Tsuruta 1964) where only the plasma 
process was taken into account as the neutrino cooling 
mechanism, we note that the faster rates of cooling by the 
URCA process and the bremsstrahlung process cause some 
important effects on the cooling behavior at higher tempera- 

g 

tures (T > 2x10 °K) , but no significant change is observed 

when T < 2x10^ °K. 
e 


OBSERVATIONAL PROBLEMS 


If a neutron star emits radiation as a black body, 
the wavelength ^ max giving the maximum intensity in the 
spectrum is given by 


t'>n\ = he/ (4.9651kT ) =0.2918/T (°K) . 

/ ) max e © 

6 7 

This simple relation indicates that when T ~ 10 to 10 °K 


the maximum comes in the soft x-ray region, 30 > \ > 3A, 

~ max ~ 

while this maximum shifts to the ultraviolet region when T 
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5 4 r. 

falls to around 10 - 10 °K» From Table 2 we see that a 

g 

neutron star with T > 10 °K has L > L (L is the lumino- 

e © © 

—5 —9 5 

sity of the sun) , while L ~ 10 L to 10 L when T ~ 10 

© © e 

4 3 5 

to 10 °Ko Also it takes about 10 to 4x10 years before a 

neutron star cools down to T ~ 10 6 °K 0 That is, a neutron 

star can last sufficiently long to allow our observation in 

the x-ray regions even though it will be too faint to be seen 

optically 0 

Until recently the observation of x-rays from outer 
space has been prevented due to the fact that the earth's 
atmosphere is strongly opaque to radiation in the x-ray 
regions. However, interstellar gases are practically trans- 
parent to x-rays of, say, X < 30l (the precise value of the 
upper limit to the wavelength depends on the distance between 
us and the x-ray emitter) „ The above considerations lead us 
to the conclusion that some of the neutron stars which are 
sufficiently close to us and which are not surrounded by 
x-ray emitting gas clouds should be observable as x-ray 
sources. 

Since it became possible to send x-ray detectors above 
the earth's atmosphere, at least 10 galactic x-ray sources 
have been discovered (Bowyer, Byram, Chubb and Friedman 1964a, 
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1964b, 1965, Giacconi, Gursky, Paolini and Rossi 1962, 1963, 
Fisher and Meyerott 1964) = Two major sources which were 
discovered first which are best known are Sco XR 1, the 
strongest x-ray source in the constellation Scorpio, and 
Tau XR 1, a somewhat weaker source near the center of the 
Crab Nebulae We use the notation of Bowyer, Byram, Chubb 
and Friedman (1965) to designate x-ray sources, A possible 
association of the strongest Scorpius source and the North 
Polar Spur, the nebular remains of a supernova explosion 
which was supposed to have occurred about 50,000 to 100,000 
years ago at a distance of about 30 parsecs away (Brown, 
Davies and Hazzard 1960), has been suggested, but otherwise 
there is no nebulosity or peculiar star in the vicinity of 
Sco XR 1. The Crab Nebula is believed to be the remnant of 
a supernova explosion which occurred in 1054 AD about 1100 
parsecs away from us. The angular size of the Sco XR 1 was 
determined to be less than 7 minutes of arc (Oda, Clark, 
Garmire, Wada, Ciacconi, Gursky and Waters 1965) ^ while the 
size of the x-ray source in the Crab Nebula was measured to 
be about 1 light year in diameter through the lunar occulta- 
tion experiment (Bowyer, Byram, Chubb and Friedman 1964b) , 
less than half the size of the optical nebula in the Crab, 
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while the optical nebula is about two to three times smaller 

-7 

than the radio size. The flux of Sco XR 1 is about 10 ergs/ 

2 -8 2 
cm -sec and that of the Tau XR 1 is about 10 ergs/cm -sec. 

The intensity of other sources is about the same as that of 

the source in the Crab Nebula . 

All of the x-ray sources, except the Scorpius source, 
lie close to the galactic plane and within 90° of the galactic 
center. This distribution resembles that of galactic novae. 
From this and other evidence it has been suggested that the 
probable origins of the x-ray sources are supernova outbursts 
(Bowyer, Byram, Chubb and Friedman 1965, Burbidge, Gould and 
Tucker 1965, Cameron 1965b, c, Morrison and Sartori 1965, 

Oda L965, Hayakawa and Matsuoka 1964, Finzi 1965). Indeed 
there is no evidence against this assumption. Many important 
supernovas have not been identified with x-ray sources of 
strength comparable with that of the Scorpius source or the 
source in the Crab Nebula. This is easily explained if 
we note that these supernovas are far more distant from 
us than the Crab Nebula or the North Polar Spur. Most 
of the discovered x-ray sources are not identified with 
known radio or optical objects. This is also no contradiction 
to the supernova hypothesis if we note the possibility 
that some of the x-ray production mech- 
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anisms will have much longer lifetime than that of the 
optical or radio emissions. 

As mentioned already, we believe that a remnant of a 

supernova explosion consists of a central condensed core in 

the form of a neutron star and of surrounding hot gas clouds 

in the form of expanding envelopes. How can we observe this 

complex assembly of matter? As far as the central core of 

neutrons is concerned there will be no hope of observation 

except as the emitter of soft x-rays in the narrow range of 

7 

about 30 to 3 a» because a hotter neutron star of T >10 °K 

e 

(corresponding to x-rays of < 3k) will cool too fast (within 
less than a day to 10 years) and a cooler neutron star of 
T g < 10 °K (corresponding to x-rays of > 30 a) will be too 
faint to be detected due to the strong interstellar absorption 
and the faintness of the star itself. The situation is more 
complicated in the expanding envelopes. 

Due to the complex nature of a supernova remnant, the 
observed x-rays can be due to non-thermal radiation from the 
hot gas clouds, the thermal radiation from the neutron star, 
or a mixture of both. The question of whether the thermal 
component can be singled out from the non-thermal component 
in the case both are present can be determined if the relative 
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strength of each component is known. In order to make some 
sensible predictions of the possible nature of some of the 
x-ray sources and their relation to neutron stars, we will 
discuss Sco XR 1 and Tau XR 1. The data on other sources are 
still too scarce for this purpose. If we know the distance d 
and the photon luminosity L of a neutron star, the flux F of 
the thermal component of x-rays reaching the region just above 
our atmosphere is found from 

(28) F = L / (4tt d 2 ) . 

If we accept the tentative association of the North Polar Spur 

and Sco XR 1, the distance and the age of the neutron star 

in Scorpius are known. Taking surface temperatures of 1 or 
6 0 

2x10 K, consistent with the soft x-ray fluxes measured by 

Friedman's group, the photon luminosity of the star can also 

be predicted o The results for our six models are shown in 

Table 11 „ The age of the star t has been taken from our 

results in the previous section (Figures 7-9) . If both the 

North Polar Spur and the x-ray source in Scorpius are indeed 

the remnants of a supernova explosion about 30 parsecs away 
4~5 

and about 10 years ago, the low mass models of M < 0.3M 

~ © 

should be excluded because they are too young. However, some 
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massive models are sufficiently old to support this hypothesis. 

The flux of x-rays from a massive neutron star is also comparable 
with the observed x-ray flux from the Scorpius source. The 
result shown in Table 11 gives rise to a possibility that the 
observed x-rays from the Scorpius source may have a predictable 
amount of thermal component, if it contains a massive neutron star of 
about 1 ~ 2 Mq. The spectral measurements of Sco XR 1 were 
reported to be inconsistent with the picture of pure black- 
body radiation. This may suggest that the nonthermal compo- 
nent from the surrounding hot gas dominates over the thermal 
component from the neutron star sufficiently to obscure the 
black-body spectrum. However, because of the relatively 
strong flux from the neutron star as calculated above, the 
thermal component in this source may be identified, if the 
remnant star is massive enough, and if longer wavelength detectors 
can be used which are more sensitive to the peak of the thermal 
spectrum. 

Next, consider the source in the Crab Nebula. Here, 

T = 910 years and we get the present temperature and luminosity 
from the results in the previous sections. The distance is 
about 1100 parsecs. From this information, the thermal 
component of the x-ray flux from a neutron star in the Crab 
Nebula can be obtained. The results are shown in Table 12, 
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for various models of different mass, nuclear potential and 
surface composition. We see that, with the exception of the 
particular model of type. Mg composition and about 2 solar 
masses, the flux is less than about 1/4 of the observed x-ray 
flux from Tau XR 1. The large size of about one light year 
already indicated that the major source of the x-rays from the 
Crab Nebula is not thermal emission from a neutron star. The 
theoretical calculation also supports this view. If a neutron 
star exists in the x-ray source, it will not be identified 
if the flux from the neutron star is too weak, as compared with 
that from the surrounding region. 


We are indebted to Mr. B. Sackaroff for assistance with 
the opacity calculations which involved the use of some parts 
of the Los Alamos opacity code in regions of temperature and 


density for which the entire code does not work. 
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TABLE CAPTIONS 


Table 1: Characteristics of our six chosen models. The models are 


identified by the type of nuclear interaction or and 

c 

their approximate mass in solar mass units. p m is the central 
density in cgs units, M/M is the mass in solar mass units, R is 


the stellar radius in km, and PC is the central pressure in 


relativistic units (PC = 


mes/cm 


6.46 x 10' 


Table 2: Photospheric properties of typical neutron stars 

(V , 0.6M ) and (V , 2M ) at different surface temperatures T . 

p O y O © 

is the photon luminosity, p^, Pp^ and x ^ are the photospheric 

o 

density, pressure and opacity, respectively, and X is the 

max 

maximum wavelength without red shift correction. 

Table 3: Atmospheric temperature distribution of neutron star 

models with M = 1M Q and R = 10 km. The temperature T at a 

specified p, and the temperature T fc and density p^ where the 

degeneracy starts are listed as a function of surface temperature 

T . 
e 

Table 4: (Models with M = 1M Q , R = 10 km) Temperature and density 

distribution near the surface at different depths (R-r) at 
given surface temperature^. The point where Ep/kT = 2.5 
indicates the thickness of the non-degenerate envelopes. 


Table 5: 


The characteristics of the hot neutron star model (V , 2M q ) 

with the Fe and Mg atmospheres. (T and T c are the surface 

and core temperatures, L , is the photon luminosity, EL is the 

P h T v 

total neutrino luminosity, a = ^ ^10 , U is the internal 

e 

energy and t is the cooling time.) 

Table 6: The characteristics of the hot neutron star model 

(V , 1.1M ) with the Fe and Mg atmospheres. (The notation is 

y 0 

that used in Table 5.) 


Table 7: The characteristics of the hot neutron star model 

(V , 0.2M ) with the Fe and Mg atmospheres. (The notation is 

y 0 

that used in Table 5 . ) 


Table 8: The characteristics of the hot neutron star model 

(V , 1.0M ) with Fe atmospheres. (The notation is that used 

0 © 

in Table 5 . ) 


Table 9: The characteristics of the hot neutron star model 


(V , 0.6M ) with Fe atmospheres. (The notation is that used 

0 © 

in Table 5 . ) 


Table 10: 

( V 


The characteristics of the hot neutron star model 
0.2M q ) with Fe atmospheres. (The notation is that used 


in Table 5 . ) 



Table 11: The observational problem of Sco XR1. (T^ is the 

surface temperature, L is the photon luminosity, t is 

ph 

the age of the neutron star of the given temperature and 
FLUX is the flux of the thermal component of the x-rays 
from a neutron star in Sco XRl, 30 parsecs away to be 
observed above the earth's atmosphere.) 

Table 12: The observational problem of Tau XRl. (T^ is the 

surface temperature of the neutron star of age 910 years, 
L ^ is the photon luminosity of the neutron star of the 
given temperature, and FLUX is the flux of the thermal 
component of x-rays from the neutron star.) 
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FIGURE CAPTIONS 


2 

Figure 1: The opacity in cm /gm as a function of density in 

3 

gm/cm , at different temperatures, calculated by the use 
of Cox's opacity code and an ionization code. Solid curves 
are for iron and dashed curves are for magnesium. 

Figure 2: Central temperature as a function of surface tempera- 

ture for models of nuclear potential V with different mass 

Y 

values. Solid curves are for iron and dashed curves are 
for magnesium. 

Figure 3: Temperature distribution near the surface of a neutron 

star. The temperature is plotted against the distance from 
the photosphere in meters for a model of one solar mass, 

10 km radius, and with various different values of the 
surface temperature T . The borders between the degenerate 
layers and non-degenerate layers are shown by crosses. 

Figure 4: Density distribution near the surface of a neutron star. 

The density is plotted against the distance from the photo- 
sphere in cm for a model of one solar mass, 10 km, and 
with various different values of the surface temperature. 

Figure 5: Energy loss rates due to neutrino pair emission from 

3 

plasma in erg/cm -sec shown as a function of electron number 



-3 

density in cm at different temperatures. 

Figure 6: Neutrino and photon luminosities and total internal 

energy of a neutron star model of 1.1 solar mass and of the 

V type nuclear potential are shown as functions of the 

Y 

core temperature of the star. The neutrino luminosity 
due to the plasma process, the URCA process and the 
bremsstrahlung process are shown separately. The photon 
luminosity for the Mg atmosphere and the Fe atmosphere are also 
shown separately. 

Figure 7: Cooling curves for the six typical neutron star models 

with iron atmospheres. Surface temperatures are plotted 

as functions of time in years. Solid curves represent the 

models with the V type nuclear potential and the dashed 

P 

curves represent the models with the type potential. 
Different models are identified by their mass values and 
the type of the potential. The points where the major 
cooling mechanism shifts from the neutrino emission to the 
photon emission are indicated by the crosses. 

Figure 8: The effect of composition on cooling curves. Surface 

temperatures of models of the type potential and with 
different mass values are plotted as functions of time in 
years for iron atmospheres (solid curves) and magnesium 



atmospheres (dashed curves). 


Figure 9: Cooling curves for models with iron atmospheres in 

the most interesting region for observation are shown in 
detail. 
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